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Abstract 

We compute the gravitational quasinormal modes of the global AdS's-Schwarzschild solu- 
tion. We show how to use the holographic dual of these modes to describe a thermal plasma 
of finite extent expanding in a slightly anisotropic fashion. We compare these flows with the 
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the elliptic flow coefficient and the thermalization time. 
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1 Introduction 

When gold nuclei are collided at y^s^vw" = 200 GeV at the Relativistic Heavy Ion Collider 
(RHIC), it is plausible that a thermal state forms in which quarks and gluons are deconfined. 
The peak temperature is estimated to be roughly 300 MeV, somewhat above the temperature 
Tc at which QCD deconfinesj^ The resulting "quark-gluon plasma" (QGP) is thought to be 
strongly coupled, partly because its collective motions are well described by hydrodynamics 
with shear viscosity much less than the entropy density: ///s <^ 1. Weakly coupled plasmas 
generally have r]/s ^ 1. Intriguingly, certain black holes in string theory which describe 
thermal states of strongly coupled gauge theories have tj/s = l/Air. Hard probes of these 
black holes, in the form of classical strings trailing across their horizons, exhibit energy loss 
which has similarities, both in magnitude and angular distribution, to experimental results 
on jet-quenching. The hope naturally arises that further and closer analogies might exist 
between real-world QGP's and string theoretic black holes. Can we "simulate" a heavy-ion 
collision, from the moment of impact to the onset of hadronization, in string theory? The 
aim of this paper is to take a step in this direction by describing expanding, cooling plasmas 
of finite extent in string theory. 

Recent accounts of the discoveries at RHIC include the authoritative summaries [ll[5l[6l[7] 
and the recent review [1]. The shear viscosity of black holes and its possible relevance to 
the QGP were developed in several papers including [HI E]- The theoretical framework for 
relating black hole physics to gauge theory is the anti-de Sitter / conformal field theory 
correspondence (AdS/CFT), formulated in [lOl [11], [12] and reviewed (among other places) in 
[131 mj. One of the first hints of the gauge-string duality was the calculation of the entropy 
density of D3-branes [TJ], which feeds into the calculations of shear viscosity to entropy 
density and also shows a 25% pressure deficit that is reminiscent of lattice results for QCD 
(see for example [T6]). 

The literature on the possible connection between string theory and RHIC physics has 
grown large. The possibility of emulating a heavy ion collision in an AdS/CFT context 
was emphasized already in [TTJ [TSl [TH], and the connection with finite-sized black holes is 
already present in these papers. Bjorken flow in an infinite medium has been studied from an 
AdS/CFT perspective, for example in [201 [21]; see also [22]. Quasinormal modes have been 

^Lattice estimates for Tc tend to fall in the range 160 — 190 MeV [1|. Recent studies [21 [3] suggest that 
there is still some disagreement within the lattice community about where in this range Tc really lies. 
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considered in [231 1211 [251 126] , the first of which especially has results that are analogous in 
significant respects to ours, albeit in an infinite static background. Hard probes were studied 
in [271 [2H1 [291 [301 [SH [32] and compared with some success to the phenomenon of jet-splitting 
[331 [M]- Quarkonium systems have also been studied: see for example [281 [SS, |36l [33 [38] . 

It is important to bear in mind that most studies of the gauge-string-RHIC connection 
have been carried out using A/" = 4 super- Yang- Mills theory at strong coupling. This theory 
is different in many respects from QCD: in particular, it does not confine, it does not in- 
clude quarks in the fundamental representation of the color group, and its vacuum structure 
is completely different. The current study is particularly closely tied to the conformal in- 
variance that A/" = 4 super- Yang-Mills possesses. Conformal invariance implies an equation 
of state p = e/3, so the speed of sound is v = The average speed of sound in the 

QGP produced at RHIC is believed to be substantially lower: for example, an average value 
f = 1/3 is used in [39]. The assumption of conformal invariance assumption may be better 
justified at the highest temperatures attained in a RHIC collision, and it is possible that a 
broader conformal regime will arise in LHC heavy-ion collisions. 

The point of entry for the present paper is the observation that the holographic image on 
Minkowski space of the global AtiS's-Schwarzschild black hole is a spherical shell of plasma 
first contracting and then expandingl^ We describe this gauge theory configuration as a 
"conformal soliton fiow" because the expanding shell has exactly the same profile as the 
contracting one, as if there were no interactions. A fiow with the same stress-energy tensor 
could be arranged in a free massless theory. It is appropriate to describe the soliton as 
conformal because its existence relies almost entirely on the ability to make a conformal 
transformation from part of S*^ x R to R^'^. In S*^ x R, the conformal soliton is simply a 
plasma at rest. 

Having established this spherical "approximation" to a heavy-ion collision in section 12.51 
we proceed to study its small deformations in linear perturbation theory. These deformations 
are known as quasinormal modes (QNM's). A large literature is devoted to them, stretching 
back to the classic computations [HI [12] for the four-dimensional Schwarzschild black hole. 
The method of choice for numerically determining quasinormal modes is usually to factor 
out the asymptotic near-horizon behavior and expand the rest in a power series around the 
horizon. Then a convergence property in the far- field limit determines u. Subtle methods 
were developed in [13], following earlier work [H], to handle the convergence issues in a 

^The dual description of global ^cJSs-Schwarzschild has been considered previously in [ID] using conformal 
transformation methods similar to those described in section [2T5l We thank G. Horowitz for bringing [40] to 
our attention. 
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numerically robust manner. A more straightforward approach, which is often sufficient for 
low-lying normal modes, was first used in the GAdSBH background in [IS] and in [IB] 
to treat gravitational perturbations of global 74(i5'4-Schwarzschild. Asymptotic expressions 
based on a different method (related to that of [Hj) were obtained in [18] for global AdS^- 
Schwarzschild, and in [19] for large black holes in anti-de Sitter space; but we are unaware 
of a full treatment of the case of interest, namely gravitational perturbations of the global 
AdS's-Schwarzschild black hole (GAdSBH). 

The plan of sections 13.11 13.21 and 13.31 is to examine in detail the tensor, vector, and 
scalar quasinormal modes of the GAdSBH and extract from them the dual perturbations 
to the conformal soliton flow (H^ . In principle, this constitutes a complete analysis of 
small perturbations. Special interest attaches to the low-lying QNM's because they tend to 
dominate late-time behavior. 

Building blocks for the quantitative study of QNM's are summarized in section [2j Sec- 
tions 13.41 and 13.51 summarize the main calculations. In section H] we consider some explicit 
examples of deformed conformal soliton flows that have some features in common with 
heavy-ion collisions. In section [5] we discuss the linearized hydrodynamic approximation to 
the low-lying quasinormal modes as well as the motion of test masses in AdS^. We also 
speculate in section O about colliding black holes in AdS^. An appendix is devoted to an 
exposition of scalar, vector, and tensor spherical harmonics on 5*^. 

The reader wishing to gain a flrst impression of our results without reading through all 
the details may wish to read sections 12. ![ 12. 5[ 12.61 13.51 14.21 14.31 and 15.51 Our estimates of 
the elliptic flow coefficient and the thermalization time can be found in sections 14.21 and 14. 3^ 
respectively, and section 15.51 is an executive summary of our results. 

2 Building blocks of the calculation 
2.1 The background metric 

The global AdS's-Schwarzschild black hole, hereafter abbreviated as GAdSBH, takes the form 



ds^ = - 




(1) 
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where dCl"^ is the hne element on a unit S^: 



dtf = gijdy'dy^ = dx^ + sin^ x{dO'^ + sin^ dcp'^) . (2) 

The background ([1]) is a solution to the vacuum Einstein equations augmented by a cosmo- 
logical term: 

4 

Rab + Jj^Qab = . (3) 



These equations of motion follow from the action 

r i9i 

(4) 



5 = ^ / d^x^g 



This action captures part of the dynamics of type IIB string theory compactified on an 
of radius L with N unit of Ramond-Ramond five-form flux, where 

L'/k^ = (N/2Tif . (5) 

AdS/CFT relates the GAdSBH background to a thermal state of SU {N) J\f = 4 super- Yang- 
Mills (SYM) on X 

The location ph of the horizon of GAdSBH is the most positive root of f{p) = 0. It is 
related to po by the following two equivalent expressions: 



For future reference, let us record the standard results for the mass, entropy, and tempera- 
ture: 



The mass is interpreted as the total energy in the boundary gauge theory (modulo a temper- 
ature-independent Casimir contribution), and it may be calculated explicitly by integrating 
the energy density (Tqq^) of fHTl) over the of radius L on which the boundary gauge theory 
resides. 

We are interested in perturbations of ([1]) which solve the linearized Einstein equations 



^In this simplest example of AdS/CFT, the five-dimensional gravitational constant is expressed in terms 

no/ 

between and k^q change. 



of the ten-dimensional one as = k^q/tt'^X^. In more complicated examples, both the (jS]) and the relation 
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following from ([3]). The perturbations should be infalling at the horizon, and they should 
have the fall-off near the boundary of GAdSBH corresponding to an expectation value of the 
gauge theory stress tensor T^n rather than an alteration of the metric (see section [23] for more 
detail). If a time-dependence of the form e"*'^'^ is assumed, then Imuj < corresponds to 
stable modes, ^ It is generally understood that QNM's dominate the approach to equilibrium 
of perturbed black holes. 



2.2 Spherical harmonics 

The strategy for solving the linearized Einstein equations in the GAdSBH background ([1]) will 
be separation of variables. So one requires generalizations of the usual spherical harmonics 
X£m(^; 0) appropriate to scalars, vectors, and symmetric tensors on S^. These harmonics are 
defined by the equations 

{V^d' + kl)S = (8) 
(ViV* + Pv)Yj = V^Yj = (9) 
(V,V* + kl)T,h = V=T,h = = T^', , (10) 

where V is the standard connection on the unit S^, indices are raised and lowered using the 
metric on the unit S^. A complete and explicit listing of the relevant harmonics is given in 
Appendix lAl following [50] . 



2.3 The hnearized Einstein equations 

Decoupled forms of the linearized Einstein equations were derived in [SI], following earlier 
work (see for example [52]). We will now briefly review the results, specialized to five 
dimensions and spherical symmetry. First, split coordinates = (r, p, x,6',0) into y°' = 
(r, p) and = (x, 0, 0). Let denote the covariant derivative corresponding to the metric 
for the "orbit space," 

dsl^-fdr' + W / = 1-4 + 4- (11) 

This metric is asymptotically AdS2- The overall strategy is to express graviton perturbations 
in terms of spherical harmonics and "master fields," denoted $, which are scalars in the orbit 

^Some works on QNM's assume a time-dependence e*"* and find stabUity when Imtj > 0. The true field 
configurations are of course made real by combining complex amplitudes with their conjugates. 
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space ( ITT!) . Let e^^ be the anti-symmetric tensor for the metric (fTT]) . such that = I- For 
tensor modes, 

Sgap = = 6gai Sgij = 2p^ Ht{t, p) Tij(x, 6, 0) (12) 

Ht = (13) 

^^^a _ Yllf}\ $ = (14) 



^t(p) = 4(¥/ + % + 4-1] . (15) 



For vector modes, 



2 ^ (16) 

= -^P' Ht{t, p) V(,V,)(x, e, 0) 

F„ = + -^9,//t = \e^/3d''{p'^''^) (17) 

jj^ga _ $ = (18) 



Mp) = iA^f--^ + k^v + ^) ■ (19) 



The quantities i^r in f|T6l) and $ in f|T7|) have nothing to do with the Ht in f|T2|) and $ in 
f lT^ . For scalar modes. 



= 2p2 [^^^1,(7-, p) §(x, 0, 0) + iyT(r, P) ^ijiX, 0, 



(20) 



S^ = S^j- = + i^,,§ (21) 

H = m + Qw w = ^ m = k'i — 3 (22) 

p2 

= J— i fa + J-daHx 

ks V 

Fa/3 = fap + D^Xp + (23) 

F = HL + lHT + -{d''p)X^ 
3 P 



(24) 
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jj^ga _ $ = (25) 



Mp) = -rirr^ h [2^^ + 36(7 - ?>w)w^ + m2(7 + hw) 



+ I2mw{llw - 9)1 + (m^ + lOSmw - 540w^) (2 - 2w - /) ' 



Again, the quantities Ht in (1201) and $ in (1241) have nothing to do with the quantities of the 
same names entering into the description of tensor and vector modes. 

The quasinormal mode calculation consists of solving the master equations (ITil) . ^ 



and (125|) . These calculations are the topics of sections |3. 1113.31 To convert these calculations 
into predictions for the boundary gauge theory, we need expressions for the holographic stress 
tensor in terms of the asymptotics of the master field near the boundary of GAdSBH. Suitable 
expressions are developed in sections YTM. 13.1.21 13.2.21 and I3.3.2"l using equations (fT2|) - (l26|) 

2.4 The holographic stress tensor 

Recall now the method for extracting the boundary theory stress tensor from an asymptot- 
ically AdS^ metric. The key ingredient is the Brown- York quasi- local stress tensor_j__ which is 
defined with respect to a hypersurface S with outward normal vector field ua 



'^AB — ^AB — KHab — -yhAB + T^^ab , , 

L 2 (27) 



9ab - nAUB Kab = -h A Venn ■ 



Here Gab is the Einstein tensor of the induced metric Kab on S. To specify S, one chooses 
a smooth, positive function g(?/) of the global coordinates which has a simple zero at 
the boundary. Then let E(e) be the locus of points where q[y) = e. The outward pointing 
normal form UAdy^ is a negative multiple of dq. The metric gab on the boundary and the 
expectation value of the stress tensor (Tab) in the boundary theory are given by 



gab = hni q'^hab (Tab) = \ lim ^T^/'^ , (2^ 



where a and b are four- valued indices along the boundary or S(e) and = SirG, the five- 
dimensional gravitational constant appearing in (j4]). 



^The sign conventions are such that Einstein's equations would take the form Gab = SttGTab if they 
appUcd (which they don't since gravity is non-dynamical in the boundary theory). Thus, Gab differs by a 
sign from the one in |53j . 
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By choosing 

ql^y^ = — (29) 

a/1 + p^/L'^ cos J + I cos X 

one arranges for Qab to be the metric for Minkowski space, although represented in the 
Penrose coordinates y°' = (r, x, ^,0)- (Section 12.51 includes an explanation of the origin of 
To transform to standard polar coordinates x™ = (t, r, 6, 0), one uses 



(30) 



L cos J + cos X L cos ^ + cos x 
The stress tensor transforms simply when passing from y°- coordinates to x"^ coordinates: 

(Tmn) is the "answer" that we want to calculate starting from a perturbed GAdSBH back- 
ground: for instance, (Tqo) is the energy density in the boundary gauge theory, and other 
components relate to the pressure, the energy flow, and so forth. 

In practice, it is onerous to work through fl27|) - fl3T]) for perturbations of the GAdSBH 
background. A more efficient route is to replace ( l29i) by q{y) = L/p, but there is a minor 
complication, discussed in detail in [53]: (Tab) as obtained via (l28!l . using q{y), now includes 
a term which accounts for the Casimir energy of the gauge theory on S''^ x R with metric 

dS^ = ~gabdy''dy^ = -dr^ + L^dQ^ , (32) 



where dCl'^ is the metric on the unit 5*^, as in ([2]). This Casimir term may be computed in a 
pure AdS^ background and subtracted away: it is independent of the black hole ^ Let the 
result of this subtraction be denoted {Tab)suh, and let 

W=\im^ = ^ = 7Ti?V4t2L2 + (L2 + r2-t2)2_ .33) 

p^oo q[y) COS ^ + COS X 

Then, as a more efficient alternative to (l3Ti) . we may use 

_ dy^dy^J_ ~ .34^ 
dx^ dx"' W"^ ' 



^Actually, it is sometimes convenient to make the subtraction work even harder: both the counterterms 
in (j27|) may be dropped, and the corresponding divergences cancel in the course of the subtraction. 
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It is useful to note that 

which follows from combining (130!) and (1331) . 

A natural expectation is that at large p, graviton perturbations of GAdSBH take the 
form 

p — iuJT 

SgAB = -^VABix, 0, 4>) + 0{p-') . (36) 

Because Sqab is a real field, ( 136|) cannot be quite right for generic, complex u. Instead, 
^9ab is a linear combination of terms of the form appearing on the right hand side, with the 
various coefficients arranged so that the sum is real. The stress tensor is a sum of terms: 

(^'afe)sub = {Tab^)svLh + {T^b^^) suh , (37) 

where the first term corresponds to the unperturbed GAdSBH background (see section 12.51 
for an explicit evaluation) and the second accounts for the perturbations. If we use (136!) and 
ignoring the reality issues for the present, (T^'^^)sub is a linear expression in tjab- 

As a further simplification, consider the case where SgpA = 0. This is axial gauge, and 
it is always possible to pass to it using diffeomorphism gauge invariance. Then using the 
foliator q{y) = L/p in fl27|) and fl28|) . one obtains 

(T,T^)sub = {Vat - rWa,) , (38) 

where (jah is the metric on S''^ x R given in fl32|) . The QNM contribution to the stress tensor 
(T^^^) on R3'^ can then be obtained from (El 



2.5 The conformal soliton flow 

To understand quantitatively how the GAdSBH background translates into a radial flow 
in the boundary gauge theory, the main steps are to justify the choice f l29|) of a foliating 
function for surfaces to which we then apply the holographic stress tensor construction. To 
this end, let us start by describing global AdS^ without a black hole: it is the covering space 
of the following hyperboloid in R^'^: 

-(X-^)2 - (X°)2 + {X'f + (X2)2 + (X'f + = -L" . (39) 
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The standard metric on AdSr, is the one inherited from R'*'^. The embedding ( !39|) is a 
convenient representation of AdS^ because the action of the conformal group 5*0 (4, 2) is 
obviousB The relation to the coordinates used in ([1]) (setting po = 0) is 



X-i = L^l + ^cos^ X° = L^l + ^sin^ = pQ, (40) 

where f2j are the components of a unit vector in R'^. Poincare coordinates are the {t, x, z) 
coordinates used to represent AdS^ in the following form: 

ds^ = ^i-df + dx" + dz") . (41) 
These coordinates cover only part of the hyperboloid and can be defined so that 



X' = L— X" = -\-l + 



y2 



z 2 V z 

where in the third equality, i runs from 1 to 3. 

The X*^ can be eliminated to obtain coordinate transformations from global to Poincare 
coordinates. To obtain explicit formulas it helps to express 

VLi = (sin X ^, cos x) ^ = (sin 6 cos 0, sin 6 sin (f), cos 9) . (43) 

f2 is a unit vector in R^. Then 

^ = — t = zJx + II? sin ^ a; = 2;-§- sin Y i1 . (a/^ 

L v^1 + pVL2 cos^ + fcosx V A^/ ^ ^ A [^y 

To understand (jS]), it helps to consider its action on the boundary. The boundary of global 
AdS^ is R X S''^ (the Einstein static universe), while the boundary of the "Poincare patch" 
covered by the coordinates {t,x,z) is R^'^ (Minkowski space). The p — oo limit of 
defines an embedding of Minkowski space in the Einstein static universe: 



sm J X sm X 



-Q . (45) 



L cos J + cos X L cos j- + cos x 
The Minkowski metric is related by a conformal transformation to the standard metric on 



^For a fuller account, see for example section 2.2.1 of 
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(A) (B) (C) (D) 

Figure 1: Four complementary views of the GAdSBH background ([1]). The first three 
are cartoons based on the Penrose diagram for global AdS^. The fourth is a true Penrose 
diagram of GAdSBH. See the main text for more detailed explanations. 

the Einstein static universe: 

ds^ = -dt^ + dr^ + r^dQ^ = W\-dT'^ + L'^dQ^) W = (46) 

cos J + COS X 

where dfl"^ is the metric on S"^ and dCl"^ = dx^ + sin^ X^^^ is the metric on S^, parametrized 
just as in ([2]). is the conformal factor. Note that Minkowski space covers the patch of the 
Einstein static universe where \t/L\ + |x| < vr (i.e. where \/W = cos ^ + cosx is positive). 

The coordinate transformation (1441) may be applied to the GAdSBH background ([T]). The 
result is a black hole which is closest to the boundary (in the sense of attaining the smallest 
value of z) at time t = 0. For t > 0, the black hole "falls" toward larger and larger values 
of z, and eventually it passes out of the Poincare patch and is no longer causally accessible 
to the Minkowski space patch of the boundary on which we wish to calculate (T^n)- The 
metric for this "falling" black hole, though entirely equivalent to ([1]), has a complicated form, 
and we will not reproduce it here. In order to evaluate (Tmn), what one does in principle is 
to expand this complicated metric around the pure AdS^, background (1411) for small z and 
extract the terms that fall to zero as z"^ for small z. {Tmn) is expressed as a linear combination 
of the coefficients of these terms: see for example (40) of [32]. In practice, because of the 
complicated form of the metric in Poincare coordinates, we prefer to work with the global 
coordinates until the last step, using flM]) or flMl) . We now recognize fl2U]) as z/L, which 
means that the surfaces S(e) on which the Brown- York quasi- local stress tensor is evaluated 
are, appropriately, surfaces of constant z (sometimes called horospheres) . 

Before turning to the computation of (Tmn), let us detour to an explanation of the four 
depictions of the GAdSBH background in figure [TJ 
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(A) The boundary at p = oo is shown as a cyhnder. Actually it is the Einstein static 
universe, R x S*^. The wedge shaped slice is the Poincare patch. The black hole passes 
into and then back out of the Poincare patch. The thin red line is the singularity at 
p = 0. The blue plane cuts the cylinder at x = Omod vr. 

(B) The cross section of figure (A) is shown, with the Poincare patch shaded in pink. 
Surfaces of constant z and constant t are indicated. Note that the black hole reaches 
its minimum 2; at t = 0. 

(C) A conformal transformation of the standard AdSs-Schwarzschild metric has the black 
hole wandering through AdS'^. The particular trajectory shown corresponds to a uni- 
form dilation in the boundary theory. 

(D) The Penrose diagram of the GAdSBH geometry. The region covered by the coordi- 
nates (r, p) is shaded light blue, and corresponds to the light blue regions in the other 
figures. The black hole singularity at p = is the upper red line. The right bound- 
ary corresponds to p = 00, and the left boundary is another copy of R x S*^ causally 
separated from the one at p = 00. 

Choosing the foliator q{%j) = Lj p a.?, described in the paragraph after (13T1) . and using (l28l) 
with q replaced by g, one may straightforwardly show that 

(Tff )sub = ^ diag I A, 1, sin^ X, sin^ x sm' ^} • (47) 

What makes this computation relatively straightforward is that the extrinsic curvature is 
simply 

Kab = - dphAB ■ (48) 

The formula fl48l) holds because of the vanishing of the shift function for the metric ([T]) 
expressed in ADM form with respect to the foliation by surfaces of constant p. 

As noted previously, there is a po-independent Casimir contribution to the stress tensor 
that has been subtracted away from the result quoted in poj) . P7|) is the expected result 
that the GAdSBH metric is dual to a stationary plasma on S^. Using 1^ and flMl) we 
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immediately extract the stress tensor on Minkowski space: 



\ mn I 



Pi (N 



) 



2 




(L2 + t2 + ^2) 




kiL^ + e + r^) 

L' + ^% 
r^L^ 



















V 











r^L^ sin^ d ) 



(49) 



Note that we are not considering in fl49l) the generahzation where an 5*0(4, 2) rotation corre- 
sponding to a dilation on the boundary is applied to the asymptotically AdS^ spacetime, as 
in figure [T]3. Such a dilation can be applied after the final answer (e.g. (149!) or a perturbation 
of it) is obtained, and it amounts, roughly speaking, to changing L in fl49p . 

As a check, we have also obtained ( H9i) more directly using the foliator q[y) = z/L and 
the formula fl3T]) . 

Some further observations about this solution may be of interest: 

1. There is no shear in the fiow fH9|) . 

2. The same solution exists in free conformal field theories, such as U{\) gauge theory. 

3. All the properties of the fiow are due to its conformal relation to a configuration on 
R X S*^ with constant energy density and pressure. So its description as the "conformal 
soliton fiow" seems appropriate. 

2.6 Estimating the size of the black hole 

Aside from new solutions generated from old ones through an 5*0(4, 2) transformation, there 
is precisely a one-parameter family of "conformal solitons" of the form (1491) : the parameter 
is po- It is more interesting to cast this parameter in a form that is manifestly invariant 
under dilations. Consider therefore the energy profile at time t = 0: 




(50) 



It is easily seen that the total energy is 



(51) 



tot — 



4L3 
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and that the average of over the probabihty distribution e(0, r) / E'tot is L^. A dimensionless 
quantity — and therefore one invariant under dilations — is the combination 

S = i?tot^ = ^||. (52) 

In our development through (jUj), L is not an arbitrary length, but instead related to the 
five-dimensional gravitational constant through ([5]). But in light of the dilation symmetry, 
we can regard L as arbitrary in flS^ . 

It is easy to argue that S should be big in order to make a comparison to heavy ion 
collisions. For example, without thinking too hard, we might take the following values for a 
central collision at RHIC: 



L = 7fm Etot = 28TeV so S ^ 10^ . (53) 

Using (152|) with A^^ = 8 (the number of different gluons when the gauge group is SU 
one finds 

There are some reasons to think that flS51) represents an over-estimate of the most appropriate 
values of S. The initial state in a Landau model of the collision is a Lorentz-fiattened pancake 
whose long axes are the extent of the gold ions (so, a radius of about 7fm, which motivated 
the choice of L in (1531) ) and whose thickness is 100 times smaller. So should we decrease 
H by a factor of 100? Maybe not by so much, because the initial state is not thermalized. 
If it thermalizes at t ~ Ifm/c, then the aspect ratio of the region containing the QGP is 
more like 10 than 100, and perhaps H ~ 10^ is appropriate. But we should also consider 
whether the value of -Etot in dSSD is appropriate. This figure comes from an estimate that 
28 TeV out of the initial 39 TeV in a gold-gold collision goes into heating and collective 
motion. But it is not at all clear that all this energy thermalizes. A conservative estimate 
of the thermalized energy would be to add up (as scalars) the transverse momenta of all the 
observed particles. This is only about 4 TeV. In quoting these numbers we have relied on 
the summary discussions in [T], which includes references to the primary literature. 

The upshot of the considerations of the previous paragraph is that S is between 10^ and 
10^, with the precise value depending on some assumptions. The corresponding range of 
Ph/ which is the parameter that will enter most frequently into QNM calculations, is 6 to 
20. Let us understand in another way what this parameter means. The peak energy density 
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in the conformal soliton flow is 



(55) 



peak 



Using the standard relations [15] for a plasma of infinite extent, 



3^8 



(56) 



one extracts a peak temperature Tpcak that satisfies 



peak 



IT \ L TT L 



2 [po 2 Ph 




(57) 



Thus Ph/L ^ 1 is equivalent to Tpcak ^ 1, which seems like a good assumption if the 
QGP thermalizes, since it is roughly the same as the condition that the Debye length should 
be much less than the extent of the plasma. A good feature of (1571) is that the dependence 
on cancels out. In fact, a overall factor could be added to the relation ([5]) between k^/L^ 
and A^, and ( 1571) would not be affected. This is because ([5]) is used in the same way to get 
f l55|) and fl56|) . In other words, fl57|) doesn't depend on the heat capacity of the boundary 
gauge theory, which is good because the heat capacity for A/" = 4 super- Yang-Mills theory 
with gauge group SU (3) is about three times bigger than for QCD@ 

3 Quasinormal modes of global Ati/Ss-Schwarzschild 

As we have remarked, the simplest gauge theory interpretation of the GAdSBH background is 
a static thermal plasma on S'^ x R. Small fluctuations in the plasma damp out exponentially 
in global time r, and quasinormal modes provide a dual description of such behavior. Readers 
familiar with [22] will note significant similarities to the results presented in this section. 
This is no accident: the in which the nearly static plasma resides is much larger than 
the plasma's reciprocal temperature, so replacing it by R'^ is at least qualitatively a good 
approximation. However, for understanding finite-sized, expanding plasmas it is crucial to 
take the additional step of conformally transforming from S*^ x R to Minkowski space. This 
introduces additional time- dependence; or, more precisely, time-dependence is expressed in 
terms of the usual Minkowski time t with respect to which the unperturbed background is 

*The heat capacity is proportional to entering into the Stefan-Bohzmann relation esB = 9*^T^- In 
the weak coupling limit, = 120 for SU{3) Af = 4, compared to 5, = 37 for QCD with two flavors and 
5, — 47.5 for QCD with three flavors. Lattice results show w 33 for QCD with 2+1 flavors at l.STc, 
whereas in the limit of strong coupling, = 90 in SU{3) A/" = 4 if one uses the standard 3/4 factor. 
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not stationary. 



3.1 Tensor modes 

In section [3. 1.11 we explain how to numerically extract the frequencies of quasinormal modes 
from the master equation (fT4|) . In section [3.1.21 we show how tensor QNM's are translated 
into traceless, conserved perturbations of the conformal soliton flow (149|) . 

There is significant overlap among the treatments of tensor, vector, and scalar modes. 
In most ways the tensor modes are the simplest. Common technical aspects are discussed 
more fully here than in sections 13.21 and 13.31 

3.1.1 Quasinormal frequencies 

Using / ^ (? and Vt ~ 15p^/4L^ for large p, and assuming e"*"^"^ dependence, one finds 
the approximate form of the master equation f[T^ : 

^1 + ^^PP'^P - ^1 <^far = , (58) 



I" 4:L\ 

which is solved by Bessel functions. The allowed solution is 



^-far = e-*-^— ^J2(^'^/p) ~ ^gg^ 

IJ UJ y/ P 



for large p. 

It is standard to introduce a "tortoise" coordinate 



dp v/I^T^tan-i + tanh-^ f 



f{p) l + 2pl/L^ 

which brings the master equation into the form 



(60) 



92 



1^ - Mp) 



$ = 0. (61) 



Because Vt(p) as p pi^, one finds immediately that the allowed solution is 

$near = e-^-(^+'^*) . (62) 

$ = e~«'^(T-''*) ig aigo a solution, but fl62p is selected by requiring that $ near the horizon is 
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a function only of r + r*, corresponding to a wave falling into the horizon. 

In the spirit of factoring out the near-horizon behavior , let us introduce a new radial 
variable, a rescaled frequency, and a separated ansatz for $: 

2/ = l-y ^ = u:L $ = e-^-(^+^* • (63) 

The differential equation for il){y) takes the form 

[s{y)dl + t{y)dy + u{y)\i^{y) = Q (64) 



:^-yff§y [{l-yff]-2iupn{l-yff 



(65) 



where 

s{y)=K{y){l-yff 
t{y) = K{y) 

u{y) = -K{y)pi,VT. 

Here K{y) is chosen so that s{y), t{y), and u{y) are polynomials in y that do not all vanish 
at any given y. The function s{y) has a simple root at the horizon, y = 0. Because y and 
kr are dimensionless, it must be possible to choose K{y) so that s, t, and u depend on pn 
and L only through the dimensionless combination pn/L. We may expand 



oo 



iy) = ^ Sii/' t{y) = ^ til/' ^(y) = ^ Mil/' , (66) 



i=l i=0 i=0 



where of course all but finitely many of t„,, and u„ vanish. A formal series expansion of 
the solution tp{y) that approaches 1 at the horizon may be developed as follows: 



(67) 



j(j - l)si+jt( 



i=0 



The strategy to determine Q is to demand that a series solution for tp{y) around y = 
should converge at y = 1. This is valid provided there are no singularities of ipiy) closer to 
y = than y = I, which is true if s{y) has no zeroes closer to y = than at y = 1. The 
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zeroes of s{y) are easily found: up to an arbitrary normalization, 

Expressions for t{y) and u{y) are more complicated, and we do not gain much intuition from 
examining their explicit forms. It is worth noting, however, that s{y), t{y), and u{y) are all 
real when Q is purely imaginary. Because of this fact, —Q* is a quasinormal frequency if Q 
is. The boundary condition $ ~ g^i^r/L p-5/2 large p implies ^{y) ~ (1 — y)^^^ for y close 
to 1, which is in contrast to the other possible boundary condition for metric perturbations, 
ip{y) ~ (1 — |/)~^/^. So if we use ip{y,^l) to denote the solution of (1671) for a specified f2, 
the boundary condition becomes = 0. This transcendental equation determines the 

quasinormal frequencies Q. 

To actually determine approximate numerical values for quasinormal frequencies, one may 
simply truncate the series solution fl67|) at some number M of terms. The resulting function 
ipMiy,^) is a polynomial in y and a rational expression in Q. So solving ipMi^,^) = 
amounts to finding the roots of a polynomial in Q. Unfortunately, for even moderate values 
of M (such as 10), the polynomials in question, when expressed analytically in terms of 
Ph/L, k^, and Q, become large and cumbersome. We have therefore found it more efficient 
to numerically evaluate iPm{^,^) on a grid of points in the complex Q plane and use these 
values to scan for zeroes. The results of the scans we performed in this way are summarized 
in tabled! Note that the dependence on the SO (A) quantum number n is fairly weak, as 
is the dependence on pn/L. For each complex quasinormal frequency Q, there is of course 
another, —Q*, as required by the reality of the master equation. 

3.1.2 Tensor perturbations of the conformal soliton 

The aim of the present section is to show that for any value of Q, a tensor quasinormal mode 
translates via AdS/CFT to a traceless, conserved perturbation of the stress tensor (T^^) of 
the conformal soliton flow (l49l) . This means that there is a priori a whole function's worth 
of freedom in deforming the conformal soliton in the direction of a given tensor spherical 
harmonic. The remarkable claim that can made using AdS/CFT is that a very restricted 
set of functions — namely, those that can be expressed as a sum of terms of the form e~*^^/^ 
where Q is one of the values found in the previous section — describes the possible behaviors 
of conformal solitons in A/" = 4 super- Yang-Mills. 

The strategy is to use (155]) to extract the contribution (T^^^)sub of the tensor QNM to 
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freq\n 


3 


4 


5 


6 


Vti/pH 


3.220 - 2.725i 
3.141 - 2.742i 
3.129 - 2.745Z 


3.264 - 2.712i 
3.151 - 2. 739i 
3.133 - 2.743i 


3.319 - 2.696i 
3.163 - 2.735i 
3.138 - 2.742i 


3.385 - 2.677i 
3.178 - 2. 73H 
3.145 - 2.740Z 




5.281 - 4.748Z 
5.193 - 4.760Z 
5.180 -4.762Z 


5.312 - 4.739i 
5.200 - 4.758Z 
5.183 -4.76H 


5.352 - 4.728Z 
5.209 - 4.756i 
5.186 - 4.760Z 


5.400 - 4.714z 
5.220 - 4.753Z 
5.191 -4.7592 




7.318 - 6.757Z 
7.216 - 6.767Z 
7.200 - 6.768i 


7.344 - 6.750i 
7.221 - 6.765i 
7.202 - 6.768i 


7.376 - 6.741Z 
7.229 - 6.763Z 
7.205 - 6.7672 


7.415 - 6.73O2 
7.237- 6. 761z 
7.209 - 6.766i 


Q4/PH 


9.350 - 8.761Z 
9.230 - 8.770Z 
9.211 - 8.771Z 


9.372 - 8.755« 
9.235 - 8.769i 
9.213-8.7712 


9.399 - 8.748Z 
9.241 - 8.767z 
9.216 - 8.77O2 


9.433 - 8.738Z 
9.248 - 8.765Z 
9.219 - 8.769Z 



Table 1: Frequencies of tensor quasinormal modes in units where L = 1. In each cell of the 
table, the top value is for pn = 6, the middle one is for ph = 13, and the bottom one is for 
PH = 20. 

the stress tensor on S''^ x R with the natural metric Qab described in f l32p . and then apply the 
conformal transformation to recover (T^^^) in Minkowski space. The full stress tensor 
for the perturbed flow is 

{Tmn) = (T^t) + (T^D . (69) 

It is from {Tmn) that we would read off quantities hke energy density that are supposed to 
be compared with the corresponding quantities in real-world quark-gluon plasmas. 
Starting from f|T2|) and f|T3|) with 

$ = e-'^r/Lp-5/2 ^ (70) 

one find that the metric perturbations can be expressed as in fj36|) with the large p asymp- 
totics 

VafS = ilai = Vij = 2Tij . (71) 

Explicit forms for Tij are given in fll48p and (11511) . The perturbations flTTl) are automatically 
in axial gauge, rjap = 0. The trace of rjab vanishes: 

1 n -iUr/L 

r'vat = -i^ff^m. = ^^^T^j = . (72) 



(Recall that gij is given by ([2]) while (jab is given by fl32|) .) So we may express fl38|) more 
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simply as 

To verify that (T^^) is conserved, first note that the connection on the factorized geometry 
X R is Vt = dr and Vj = Vj, where as before Vj is the standard connection on the unit 
S^. Thus 

^a^^QNM^ = -dAfT"") + ;^V^(f;r'') • (74) 

The first term vanishes because {T^f^^) = 0. The second term is proportional to V*Tjj, 
which also vanishes according to the definition ( |TOl) . 

Putting flM|l . fl7T|) . and fl73l) together, one winds up with the following final formula for 
a tensor perturbation to the conformal soliton flow: 

where Q is any of the dimensionless complex quasi- normal frequencies found in section 13.1. ![ 
and rjab is given as in fl7Tl) PI Note that, after plugging in formulas for T.^ from (I148P or (11511) 
and using (1551) and (15S]) . the components of (T^^^) arising from fl75]) are explicit alge- 
braic expressions in terms of t, r, and trigonometric functions of 6 and 0. Unfortunately 
these expressions are fairly complicated even for the simplest perturbations, and we will not 
reproduce them here. 

The alert reader may notice that what we demonstrated in (TMll is not the conservation 
of (T^^^) with respect to the standard, flat connection on Minkowski space, but instead 
the conservation of {Tab)suh with respect to the Christoffel connection Va on x R, with 
metric as given by fl32|) . But, together with tracelessness, conservation carries over from one 
conformal frame to another: one may think of the power of W in fl75]) as just the right factor 
to ensure that the conformally transformed stress tensor is conserved with respect to the 
Christoffel connection of the conformally transformed metric. 

3.2 Vector modes 

Vector quasinormal modes are qualitatively similar to tensor modes except that, in addition 
to a tower of modes that is nearly isospectral with the tensor quasinormal frequencies shown 
in table m, there is a single low- lying mode whose rescaled frequency Q = ujL scales as L/pn, 
rather than the usual pn/L. Discussion of technical aspects in common with the case of 

minor subtlety arises for the case where ^ is purely imaginary, to be discussed in section [?.2. 21 
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tensor QNM's will be abbreviated. 



3.2.1 Quasinormal frequencies 

Using / ^ / L'^ and Vy ~ 3p^/4L^ for large p, and assuming dependence, one 

finds two asymptotic solutions to the master equation f|T8l) . with p^/^ and p~^/^ behavior, 
respectively. The former corresponds to metric deformations, which we disallow. The latter 
may be developed as a power series: 



~8 ^ 



0{p 



(76) 



Expansion of $ near the horizon in a power series times an overall exponential proceeds in 
a fashion identical to the treatment (!58|) - (l67|) for tensor modes, except that Vt is replaced 
by Vy, which also vanishes at the horizon. Even the functions K{y), s{y) and t{y) are the 
same, though u{y) and hence ip{y) are different. Thus the series expansion method described 
below fl67j) is justified, and we used it to scan for quasinormal frequencies Q = ujL in the 
complex plane. 

A notable feature of the vector quasinormal frequencies is the occurrence of a single 
low-lying mode for each allowed value of fcy, namely ky = ti^ — 2 for n = 3, 4, 5, . . ., with 
frequencies which appear to be a good fit to the following formula, similar to a formula found 
in [16] in the AdS/^ case: 

^^-\-{n'-A). (77) 
4pH 

(See section [3^ for a more systematic discussion of the best fit formulas to low-lying quasi- 
normal modes.) Actually, each VL in ( iTTl) corresponds to a pair of {ri^ — l)-dimensional 
irreducible representations of 0(4), filled out by choosing to be one of the explicit forms 
given in (fT43|) or ffT46|) . 

Continuing the formula (177|) to n = 2 would suggest that there are vector zero modes, 
forming an adjoint representation of 0(4). This is obviously true: there are rotational moduli 
of the black hole, and they indeed fall in the adjoint of 0(4). But, as the authors of [51] 
note, the derivation of the master equation f|T6|l relies upon assuming n > 3. For the special 
case n = 2, a simpler treatment is possible and was given in [52l [5T] . The relevant Einstein 



^"^As with tensor modes, (|18p can be solved in terms of hypergeometric functions in the case of pure AdS^,, 
but wc will not need this solution. 
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equations may be cast as conservation of a "field strength:" 

d, (p^e"^F„^) = wfiere F^p = pdj-^ -pdp^. (78) 

Fai3 has nothing to do with the symmetric tensor Faf3 entering into (12^ . There is only one 
solution of fl75]) consistent with axial gauge, /p = 0, and the vanishing of fr for p ^ that 
is imposed on us by demanding that the metric on the boundary is not perturbed. It is 



where Q is the angular momentum resulting from turning on the perturbation. Unlike the 
n > 2 cases, there is no tower of non-zero quasi-normal frequencies. The zero-mode (1791) is 
the only allowed perturbation in the n = 2 sector. Setting Q = 1, one immediately obtains 

Vrr = Vri = \^i V^J = ^ ■ (80) 

Besides the low- lying mode described in ( !77|1 . there is a main series of quasinormal fre- 
quencies for all modes with n > 2. Remarkably, these are nearly isospectral with the tensor 
modes in the cases that we explored numerically, namely pn = 6, 13, and 20 with 3 < n < 6: 
see table |2l 

3.2.2 Vector perturbations of the conformal soliton 



In this section, we wish to show that the expression f[73l) still applies, with suitably altered 
rjab, and that the resulting stress-energy tensor is conserved. 
Plugging (176!) into f|T7|) we get 

Fr = ^^l^Lt^e-f^-Al + O (p-^) F, = -^fie"^^^/^4 + O (p~') . (81) 
4 p-^ ^ p^ ^ ' 

The quantities and Ht are not separately gauge invariant. Roughly speaking, gauge 
choice at the level we are interested in amounts to choosing Ht- Axial gauge, fp = 0, can 
be arranged by choosing 

4 p"^ ^ 

Then 

/. = ^e-^-/^l + 0(p-^) (83) 
4 p-^ 
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freq\n 


2 


3 


4 


5 


6 


Ph^o 





-1.255i 
-1.246i 
-1.237z 


-3.038i 
-3.002i 
-2.990Z 


-5.369i 
-5.269i 
-5.247z 


-8.283i 
-8.053i 
-8.01H 


Qi/Ph 


N/A 


3.184- 2. 736i 
3.133 - 2.744i 
3.125 - 2.746Z 


3.207- 2.729i 
3.138 - 2.743Z 
3.128 - 2.745Z 


3.235 - 2.720Z 
3.145 - 2.741Z 
3.130 - 2.744Z 


3.268 - 2.708Z 
3.152 - 2.738Z 
3.133 - 2.743Z 




N/A 


5.256 - 4.755i 
5.188 -4.761i 
5.178 -4.762i 


5.272 - 4.750Z 
5.192 -4.760Z 
5.179 -4.762i 


5.292 - 4.744z 
5.196 -4.759Z 
5.181 -4.76H 


5.316 - 4.737i 
5.202 - 4.758i 
5.183 -4.761i 


^■i/pH 


N/A 


7.299 - 6.76H 
7.212 - 6.767i 
7.199 - 6.768Z 


7.312 - 6.758Z 
7.215 - 6.766Z 
7.200 - 6.767Z 


7.328 - 6.753Z 
7.219 - 6.765Z 
7.201 - 6.767Z 


7.348 - 6.747Z 
7.223 - 6.764i 
7.203 - 6.767Z 


VLa/Ph 


N/A 


9.334 - 8.764Z 
9.227- 8. 769z 
9.211 - 8.7702 


9.345 - 8.761Z 
9.230 - 8.769i 
9.212 - 8.770i 


9.359 - 8.757Z 
9.233 - 8.768Z 
9.213 - 8.769i 


9.376 - 8.752i 
9.237 - 8.767Z 
9.215 - 8.769i 



Table 2: Frequencies of vector quasinormal modes in units where L = 1. In each cell of the 
table, the top value is for pu = 6, the middle one is for pn = 13, and the bottom one is for 
Ph = 20. The low-lying frequencies, approximated by (!77|) . are shown as VLq. The n = 2 
column contains only the zero mode described in (!78!) - (l80l) . 

and 

2 — k"^ LVt 

Tlrr = ^ r/r^ = ^^Vi ^/i,- = -Z — V(iV,) . (84) 

The trace g"'''riab is proportional to ^*-' V(jVj), which vanishes by the definition ([9]) of Vj. So 
we can again use the simplified relation ( 1731) for the stress tensor, and what remains is to 
check that e^'^^'^^^rjab is conserved. For h = t this is easy: 

and the right hand side vanishes because it is proportional to V*Vj. For b t we have 

2-en n ^^^^ 

4 L 4L ^ ' 
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To show that the last expression vanishes, we compute 



(87) 

-klNi - Rh'^j^k = -klNi + Ri^Nk = (2 - kl)Ni 



1 1 



where in the last step we used Rij = 2gij for the unit S^. 



3.3 Scalar modes 

Scalar modes are for several reasons more technically challenging than the tensor and vector 
cases: 

1. The potential Vs{p) from (126|) is more comphcated, resulting in polynomials s{y), t{y), 
and u{y) that are higher order than before. 

2. The far field behavior of tp{y) (defined analogously to fl63|) ) is either y^l — y or ^1 — y x 
log(l — y).The latter is disallowed by the usual conditions of not changing the metric on 
the boundary, but it is numerically challenging to tell the difference between functions 
with these two behaviors. In particular, we can no longer simply evaluate ipNiu) at 
y = 1 and look for the values of Q where it vanishes: with sufficient accuracy, it will 
vanish for all Q\ 

3. There are two nearly degenerate quasinormal frequencies close to the origin, and it is 
a challenge to obtain enough numerical resolving power to distinguish them. 

As a final illustration of the challenges, we note that the methods of sections 13.21 and 13. 
applied without regard to point 2 and carried to 120 terms in the series and 150 digits 
of working precision, sometimes gave values of Q with Imf2 significantly greater than 0, 
contracting the stability proof of [5^ . 

To surmount the difficulties just described, we adopted a modified approach illustrated 
in figure [2l The plan is to use numerical integration of the differential equation fl64|) (im- 
plemented through Mathematica's NDSolve) to compute ip{y) from initial conditions at 
y = yi = 1/4 obtained through a power series expansion of the form fl671) around y = 
up to a value y = yj that is half as far from y = 1 as the zeroes y = yi and y = yi. In place 
of we computed the Wronskian between the numerically evaluated ip{y) and the series 
expansion of the allowed solution around y = 1. This Wronskian vanishes only if the same 
il^iij) that has regular behavior at the horizon has ^J1 — y behavior near y = 1: in other 
words, for quasinormal modes obeying standard AdS/CFT boundary conditions. 
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A 



y 



(ii) (iii) (iv) 

12 



®y2 

Figure 2: The zeroes of s{y) are shown as crosses in the complex y plane. The red cross 
at the origin is a single zero; all others are double zeroes. The region (i) is where a power 
series expansion around y = was used in computing 4j{y). The line (ii) is where ip{y) 
was computed by numerical integration, seeded with initial conditions at the boundary with 
region (i). The region (iii) is where ip{y) was computed using a power series around y = 1. 
The powers in this latter series are half-integral, so there is a branch cut which may be 
located as shown by (iv). 
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The rationale for this plan is that ip{y) must be analytic away from the zeroes of s{y) and 
(sometimes) branch cuts emanating from these zeroes. Up to an arbitrary normalization, 



s{y) = y{y - lf{y - 2f{y - y,f{y - y{)\y - y^fiy - mf 




(88) 



Further details of the computation of quasinormal modes are presented in section 13.3. 1[ fol- 
lowed by a discussion in section [3.3.2l of the dual stress tensor perturbations of the conformal 
soliton flow. 

3.3.1 Quasinormal frequencies 

From the large p asymptotics / ~ P^/L'^^ Vs ~ 1^^^-, ^-nd the assumption $ oc e"'*^"^/-^, 
one finds that the two linearly independent solutions of fl25l) behave asymptotically as 
p^^/^iQgp and at large p, corresponding to the \/\ — y log(l — y) and \/\ — y be- 

haviors for ^/'(y) noted previously. The former solution should be disregarded because it 
corresponds to deformations of the boundary metric. In section 13.3.11 we will have occasion 
to use the series expansion 



where m = A;| — 3, as in (j22l) . 

After defining ^/'(y) as in (l63l) . a treatment entirely analogous to (1671) is possible near the 
horizon. A significant difficulty is that the recursion relations involve many terms, so that 
evaluating aj symbolically is impractical for j greater than a few. Numerical evaluations of 
the Oj lose accuracy as j increases on account of the round-off error arising from iterated 
evaluations of the recursion relations. Experience shows that high accuracy can easily be 
obtained from the power series expansion up to y of 1/4, or even 1/2. 

The series expansion of the allowed ^\ — y solution around y = \ proceeds in a fashion 
rather similar to (1671) . and we omit the details. Experience shows that high accuracy can 
easily be obtained from this series expansion out to 1/2 or 2/3 of the radius of convergence. 
Unfortunately, as illustrated in figure [2l the two regions of good convergence do not overlap. 



iHr/L -1/2 ^ 



+ 




m(l + 
32 



(89) 
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freq\?T, 


3 


4 


5 


6 


Re Qq + ph Ini 


1.652 -0.826i 
1.637 -0.832i 
1.640 -0.842Z 


2.292 - 1.972i 
2.248 - 1.994i 
2.241 - 1.997i 


2.949 - 3.428i 
2.854 - 3.485i 
2.839 - 3.494Z 


3.633 - 5.179i 
3.462 - 5.303i 
3.435 - 5.321Z 


Qi/Ph 


3.172 - 2.740Z 
3.131 - 2.745Z 
3.124- 2. 746i 


3.188 - 2.735Z 
3.134 - 2.7442 
3.126 - 2.746Z 


3.208 - 2.729Z 
3.138 - 2.743Z 
3.127- 2. 745z 


3.233 - 2.723Z 
3.143 - 2.741Z 
3.130 - 2.744Z 




5.247 - 4.758Z 
5.186 -4.762Z 
5.177 -4.763i 


5.258 - 4.7552 
5.188 -4.762Z 
5.178 -4.7632 


5.272 -4.751Z 
5.191 -4.761Z 
5.179 -4.762i 


5.290 - 4.747z 
5.195 -4.760Z 
5.180 -4.7622 




7.291 - 6.765Z 
7.210 - 6.768Z 
7.197- 6. 769z 


7.300 - 6.762Z 
7.212 - 6.768i 
7.198 - 6.769i 


7.311 - 6. 759z 
7.214 - 6.7672 
7.199 -6. 769z 


7.325 - 6.756Z 
7.217- 6. 767z 
7.200 - 6.768Z 




9.327- 8. 768z 
9.225 - 8.772Z 
9.209 - 8.772i 


9.335 - 8.766Z 
9.226 - 8.771Z 
9.209 - 8.772i 


9.344 - 8.7632 
9.228 - 8.771Z 
9.210 -8. 772i 


9.356 - 8.760Z 
9.231 - 8.7702 
9.211 - 8.772i 



Table 3: Frequencies of scalar quasinormal modes in units where L = 1. In each cell of the 
table, the top value is for pn = 6, the middle one is for pn = 13, and the bottom one is for 
PH = 20. 



To cross the gap between these two regions, we resorted to numerical integration of the differ- 
ential equation for tp, of the form fl6^ (but, obviously, with appropriately altered polynomials 
s{y), t{y), and u{y)). As previously remarked, the function of Q whose zeroes determine the 
quasinormal frequencies is, in our approach, the Wronskian aty/ = l — |1 — ?/i|/2 between a 
numerically integrated form of ipiy) which obeys the usual horizon boundary conditions and 
the allowed series solution around y = 1. 

Similarly to the vector mode case, the scalar modes consist of low-lying modes and a 
main serieso Remarkably, the main series turns out to be again nearly isospectral with the 
tensor and vector modes for the cases that we examined, namely for pn = 6, 13, and 20, and 
3 < n < 6: see table O The low-lying scalar modes differ from the vector ones in that they 
have non-zero real parts. Their imaginary parts are roughly 2/3 times the imaginary parts 
of the corresponding vector modes. 



^^We believe that the same pattern should persist for scalar modes (more commonly described as even 
modes) in AdSi. Main series modes were found in [46], but the low-lying modes were not. In section [Ql we 
will comment on the hydrodynamical interpretation of the low-lying modes. The hydrodynamical discussion 
should apply with small adaptations to the AdS^ case. 
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3.3.2 Scalar perturbations of the conformal soliton 

Analogously to sections 13.1.21 and 13. 2. 2^ we wish to show in this section that fl73p can be 
used to compute the stress tensor dual to a scalar quasinormal mode, with suitably altered 
rjab, and that this stress tensor is conserved. 
Plugging (jSnD into fl21|) gives 

= ^ + O (90) 

2 p-^ ' 

where 

^ ^ kl{3kl - 8) - lon^kl - 2) + 7n' ^_,^^/, _ ^9^^ 

32 

We also note that (12^ implies the following relations between Fajj and F (see [51] for a more 
detailed explanation): 

° + 2F = (pF,^ - 2pF(?,^) = . (92) 

The above equations (lUO]) and (19^ . together with the axial gauge conditions fp = fpp = 
frp = 0, determine unambiguously the leading order dependence on p in the metric pertur- 
bations fl20|) . Namely, fr, frr, Hl, and Ht are: 

12 ^ 
^, = MRz^e---/^-i + 0(p-) 

48 p4 ^ 

When combined with (l36l) . the above expressions give explicit formulas for r]ab- 

kUkl-3)^ ^ zmks{kl-3)^ 



Vtt = — — s n 



12 ' 12 

2 

S 



LV,^ , 1 (94) 
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{kl-3)g,,S+^{kl-3n^)% 
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which in turn can be used to find the VEV of the QNM contribution to the stress tensor in 
the boundary theory via 

Equation fl38|) contains the trace g'^^rjab, which was zero for the tensor and vector modes. 
Likewise, for the scalar modes, we obtain from (19^ 



-aft _ 3) kl 

9 nab- §+- 



3(A;|-3)§+-(/t|-3fi2)§^ 



0. 



(95) 



We have used S\ = 0, which follows from the definitions (12T]) of Sjj and ([8]) of S. Therefore, 
?7 is traceless in this case too, and we can again use the simplified equation (173|) to compute 
the VEV of the stress tensor. 

It remains to check that the holographic stress tensor is conserved, which is the same as 
checking that V" (e-^^^/-^r/afe) = 0. To do this, we will use the identities [52] 



as = 



2 ^1 3 J 

3 ks 



(96) 



The first of these is trivial, and the second and third follow in a more or less straightforward 
way from the definitions ^ and ^ of §, §i, and S^^. To check that V''(e-^^^/-^r/a^) = 0, 
we compute as follows: 



-iQr/L 



:.^ fe|(fe|-3) g _ 1 zLnks{kl-3) ^,^ 



12 



L2 



12 



which vanishes by the second relation in (1961) . Similarly, 

V" (e-*^^/^r/„,) = -9. (e-^^^/^r/„) + -^V^' (e-*^^/^r/,,) 

_ -iflr/L 



-iflr/L 



12 ' 36 24 

l]2fc5(fc|-3) kUkl-3) , fc5(fc|-3(]2)(fc| -3) 



12 



36 



36 



(97) 



(98) 



To get to the last line (which obviously vanishes), one uses the first relation in (1961) on the 
second term of the previous line and the third relation in (!96|) on the third term. 
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3.4 Approximate formulas for the quasinormal frequencies 

We find that all the quasinormal frequencies for the main series of the tensor, vector, and 
scalar modes are approximately described by 

1] = fii + Afi(g - 1) q>l 
VLi/ph = (3.195 ± 0.009) - (2.743 ± 0.001)z (99) 
An/pH = (2.020 ± 0.003) - (2.0061 ± 0.0005)^ 

for ph = 6, 13, and 20, and 3 < < 6. (Note that we are still setting L = 1.) Here, the 
uncertainties in the regression coefficients are simply the standard errors that we obtain from 
a simultaneous fit of all the frequencies from tables [H [21 and |3l Inspection of these tables 
shows some trends not captured in (!99l) . for example weak dependences on n and phlI 

The temperature as measured with respect to global time is Tg ^ 2ph/tiL'^ for p^ ^ L. 
See d?]) for a form that is valid for arbitrary pu/ L. It is striking that /WL/2TxTgL ^ 1 — i. 
If we revert to the notation uj = Q/L, then the imaginary part of Au is approximately the 
first Matsubara frequency. This clearly calls for some interpretation in terms of a Euclidean 
continuation. 

The low-lying vector modes are approximately given by 

n = -^{n^ -nl) a = 1.018 ±0.008 = 4.14 ±0.17 (100) 

for ph = 6, 13, and 20, and 3 < ri < 6. For the same range of n and pn we find that the 
quasinormal frequencies of the low-lying scalar modes are given by 



= 1.031 ±0.012 n^^^ = 1.46 ±0.41 
= 0.985 ± 0.007 i = 3.91 ±0.14. 

An independent calculation of the low-lying vector and scalar modes can be done using 
linearized hydrodynamics on 5*^ x R, which we will do in section 15.11 Formulas fllOOl) and 
(110 II) are written in a form that makes easy the comparison to predictions from linearized 
hydrodynamics, fll27p and (11231) . 



^^An example of such a trend in the scalar modes is the slight increase in the magnitudes of Refio and 
p/f ImOo as ph increases. An apparent exception is the real part for the n = 3 case. Closer inspection of 
the Wronskian function in this case shows that the numerics arc not as reliable as usual for ph/L = 20. 
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Figure 3: The real (left) and imaginary parts (right) of the low-lying quasinormal frequencies 
for the scalar modes with 3 < n < 11 (we have set L = 1). The green diamond-shaped points 
correspond to pn = 6, the red triangular ones to pn = 13, and the blue star-shaped ones to 
pH = 20. The dashed line shows the predictions from linearized hydrodynamics, as explained 
in section [5?T1 The fit (110 II) was done for the leftmost four points shown (3 < n < 6), but 
simultaneously for all three values of pu — so 12 data points altogether. These points are the 
ones inside the inset frame in each plot. The fit (I102p comes from only the pu = 20 points. 



We extended our calculations of the frequencies of scalar QNM's to higher n and found 
that deviations from (110 II) are fairly significant, especially for pu = 6: see figure [31 On the 
other hand, for = 20, the fitting function in (1 10 II) works well for 3 < n < 11, but the 
best fit parameters change a little: 

= 1.0170 ±0.0021 1.41 ±0.16 

(102) 

tti = 0.9922 ± 0.0006 nl ^ = 3.86 ± 0.04 . 

As is evident from figure [31 the pn = 13 frequencies are fairly close to the values for p^ = 20. 

The overall pattern of quasinormal frequencies is similar to that found in [23]. As noted 
earlier, this is not accidental: large pn/L means that replacing by is for many purposes 
a good approximation. 



3.5 Summary of results 

Let us summarize in a single formula all that we have learned about the stress tensor profiles 
dual to perturbed GAdSBH's: 

(Trnn) = («) + E (^Z^)'^""^ ^"^^""^^^ ' (103) 

R 
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where Qr are complex constants. The rest of the notation in fllOSp is explained in the next 
paragraph. The expression fll03p is a large family of exact solutions to the equations 



V™(T„„) = (T™„) = 0. (104) 

These solutions are dual to approximate descriptions of deformations of GAdSBH — approximate 
because we used the linearized Einstein equations. It is worth reiterating a point from the 
end of section 13.1.21 after plugging in formulas for the spherical harmonics developed in ap- 
pendix [X] and using fl33|) and fl35l) . the components of (T™) are explicit algebraic expressions 
in terms of t, r, and trigonometric functions of 6 and 0o 

In fll03p . (T^'J) is the unperturbed conformal soliton flow (HOl) . A^^ is the dimension of the 
gauge group, L can be regarded as an adjustable parameter (because of dilation symmetry), 
and W is the conformal factor defined in (l33l) . i? is a "composite index" that labels a 
particular basis vector in a particular representation of 0(4) and a quasinormal frequency 
appropriate to that representation. For convenience we also include a choice of pn/ L in the 
information that R specifies, but actually this choice is fixed once we select the background 
flow (T^^). In other words, pn/L is the same for all the terms in (I103p . Consider for example 



R = {tensor, even, n = 5, £ = 3, m = -1, pu/L = 20, VL ^ 66.4 - 53.9z, p = 0} . 

(105) 

The first three elements of this list select the 0(4) representation as the symmetric part of 
7l ® 3r © 3l C?) 7r (see the discussion around fll42p ). The next two elements of the list 
select the basis vector corresponding to T^J™(5, 3, —1). (See fll48p for an explicit expression 
for T°J°'^(5, 3, —1).) To form 7]^^ we refer to fl7T]) . For vector or scalar modes we would refer 
instead to (l84p or (1941). 

Once the representation and pu/ L are selected, there are discretely many possible quasi- 
normal frequencies, namely the roots of iI){\^VL) where ip{y,Q) is a solution of the modified 
master equation (1641) (note however that a more sophisticated prescription is needed in the 
case of scalar QNM's, as described in section 13.31) . In sections I3.1.H 13.2.11 and 13.3.11 we 
made a numerical determination of some of these frequencies, as summarized in tables (H 
|2l and [3] and the approximate fitting forms discussed in section 13. 4[ If a particular root of 
has some multiplicity J > 1, then the allowed time dependences at this frequency 
are r^e~*^^/^ where < p < J — 1. These time dependences can be thought of as arising 

^■^An exception is the case where some p{R) is non- vanishing, in which case powers of logs of rational 
expressions in t and r would also be involved. However, as noted in the main text, non-vanishing p{R) is 
only a theoretical possibility, and we did not encounter any examples in actual computations. 
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from differences of pure exponentials witli infinitesimally separated frequencies Q. A familiar 
example is the critically damped simple harmonic oscillator. Although we have included in 
(11031) the possibility of non-zero p, our numerical explorations yielded only single roots of 

The case of vector perturbations in the adjoint of 0(4) (n = 2), as remarked around (I78|l . 
is special in that there is only one mode, and it has = and ?7m„ given by fISU]) . 

The only input from black hole physics (at the level of the current discussion) is the 
particular values of Q that we must use. The spectrum of possible Q is symmetric under 
Q — > —Q*, so by appropriately choosing the complex coefficients Qr one may ensure that 
{Tmn) is real. Although (Tmn) is exact/?/ traceless and conserved for any choice of the (and 
indeed for arbitrary Q), in general it might have pathologies like negative energy density at 
sufficiently early times (e.g. before the "collision"). Thus we should regard the perturbed 
conformal soliton as a framework best justified for addressing the late-time behavior of 
colliding conformal matter. To achieve a better description of early times via AdS/CFT, it 
would be better to avoid the linearized approximation to Einstein's equations and consider 
some physical collision in anti-de Sitter space which results in a black hole that "rings 
down" through QNM's at late times towards GAdSBH or some other stationary solution. 
For example, one could consider colliding shock waves, as advocated as a description of RHIC 
collisions in [1^, or colliding black holes. 

4 Case studies of perturbed conformal soliton flows 

The aim of this section is to take the extended technical treatment of quasinormal modes 
of the GAdSBH background ([T]) from section 13.51 and attempt to extract some flows in the 
boundary gauge theory that more closely resemble heavy-ion collisions than the original 
conformal soliton (jUj). There is a problem with this program at the level we are able to 
approach it in this paper: it's not clear a priori that the late-time regime in which QNM's are 
a good description of the black hole dynamics overlaps with the early time regime in which 
the temperature is high enough for conformal invariance to be an approximate symmetry of 
the QGP. The situation might be improved somewhat by going to higher temperatures on 
the experimental side, as LHC experiments will do. On the theoretical side, it is possible 
in principle to go beyond the QNM approximation with numerical simulations of the full 
Einstein equations. It is worth recalling that although a real-world QGP at RHIC would be 
fully hadronized by times t ^ L/c = 7fm/c, a truly conformal plasma (such as a A/" = 4 
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super- Yang-Mills plasma) would simply continue to cool and expand without going through 
any phase transitions. 

Given the limitations just described, it is best to think of the studies below as preliminary 
forays intended to help define the questions we might ask of more substantial explorations. 



4.1 Slow vector modes 

Consider the low-lying vector modes whose quasinormal frequencies are given approximately 
in fl77p . together with the zero modes described in the subsequent paragraph. Simply for the 
purpose of getting started, let us consider the mode specified by 

Ri = jvector, even, n = 3, i = 2, m = 0, ph/L = 20, p = o| . (106) 

From now on, we will mostly set L = 1, which is equivalent to measuring all lengths in units 
of L = 7 fm. The energy density and the radial Poynting vector are 



e = (Too) = e'^^t, r) + Qi en, (t, r, ^, 0) + . . . 
Sr = -{Tor) = 5f^t,r) + Qi SrMt,r,e,<p) + 



(107) 



Choosing = 8 X 10^ and A^^ = 8 results in a flow illustrated in figure ffl^j This flow is 



reminiscent of a Lorentz-fiattened pancake expanding in something resembling Bjorken flow. 
Alternatively, if the z axis is replaced by the y axis so that the views in figure H] are along the 
beamline, the flow appears somewhat similar to elliptic flow, where an initially anisotropic 
state expands faster in along its shorter axis. 

The Bjorken flow analogy of the previous paragraph might contain a grain of truth, but 
we believe that the physics of the solution (I107P is wrong for describing elliptic flow. Consider 
the following points regarding the solution (11071) : 

1. The energy density is spherically symmetric at t = 0. This feature is common to all 
flows constructed by deforming the radial conformal soliton only by vector and tensor 
perturbations. 

2. The radial motion is, on average, inward rather than outward between t = —0.5 and 
t = 0. After t = the radial motion is mostly outward. (Recall that we have set L = 1. 
If we restored dimensions with L = 7fm, t = —0.5 would translate into t = — 3.5fm/c.) 



^^Choosing Qi = 10^ A'^^ for any N results in precisely the same flow up to an overall normalization. 
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Figure 4: Snapshots of the energy density of the flow specified by (11071) with Qi = 8 x lO'^. 
The z coordinate measures position along the beamline (where Q = or tt), while the x 
coordinate is transverse. The flow is azimuthally symmetric, which is to say symmetric 
around the z axis. Note that the scales of the x and z axes change from frame to frame. 
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3. There is negative energy density for times earlier than approximately t = —0.56. If 
Qi were increased, there would be negative energy up to a later time. Negative energy 
density at sufficiently early times is common to most linear perturbations fll03p of 
the conformal soliton flow, because the perturbations are signed quantities that grow 
exponentially as t ^ — oo. This pathology is presumably cured when one goes beyond 
linear perturbation theory in the GAdSBH background. 

4. Referring to the t = —0.5 snapshot, one sees that the energy density profile is flattened 
in the z direction by at best a factor of a few. Real gold-gold collisions start with an 
aspect ratio possibly as high as 100, although by the earliest time that thermalization 
is expected, this aspect ratio is probably closer to 10. 

From the earliest time {t ^ —0.56) that the linear approximate can be used without en- 
countering negative energy density, the radial component of the Poynting vector is more 
positive on the horizontal axis of figure H] than on the vertical axis. Moreover, Q is purely 
imaginary. So instead of describing some compressed state accelerating preferentially along 
this horizontal axis, it seems that we have described a flow that starts with outward veloci- 
ties preferentially horizontal and that manages to keep expanding horizontally despite some 
modest damping. This is not how elliptic flow is thought to arise. 



4.2 Slow scalar modes and elliptic flow 

Now consider the low-lying scalar modes whose quasinormal frequencies are given approxi- 
mately in fllOip . To get started, let us consider the mode specified by 

i?2 = {scalar, n = 3, i = 2, m = 0, pn/L = 13, ^ 1.637 - 0.064z, p = 0} . (108) 

We will again set L = 1 in the rest of the discussion. In analogy to (11071) . we exhibit an 
explicit special case of (11031) as 

e = (Too) = e^^^(t, r) + 2 Re {Q2 en, {t, r,9,(P)} + ... 

(^109) 

Sr = -{Tor) = ^t, r) + 2 Re {Q2 Sr^t, r, ^, 0)} + . . . . 

We must explicitly take the real parts in fll09p because the quasinormal frequencies are 
complex. This is in contrast to (I107p . where the frequencies are purely imaginary and the 
quasinormal mode expressions are already real. Choosing Q2 = —10^ and A^^ = 8 results 
in a flow illustrated in figure O (As before, the same flow up to a normalization would 
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Figure 5: Snapshots of the energy density of the flow specified by fllOOp with Q2 = —10''. 
The z coordinate measures position along the beamline (where 6 = or tt), while the x 
coordinate is transverse. The flow is azimuthally symmetric, which is to say symmetric 
around the z axis. Note that the scales of the x and z axes change from frame to frame. 

result from choosing any Q2 and A^^ with Q2/N'^ = —10^/8.) Once again, we can attempt 
an analogy to the longitudinal collective motion of a central collision, or, by changing axes 
z y so that the plots in figure [5] are views down the beampipe, elliptic flow in a non-central 
collision. Both analogies may have some merit, but let us focus on the possible connection 
to elliptic flow. At t = 0, the radial and tangential components of the Poynting vector are 
small: typically on the order of 10% or less of the change \Q2^R2\ in the energy density. 
This is in stark contrast with the situation for vector modes, where e^^ = at t = while 
the Poynting vector is non-zero. The smallness of the Poynting vector means that we are 
starting at t = with an approximately stationary plasma that has been compressed into 
an elliptical shape. Its preferential expansion along the short axis of the ellipse must then 
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be due to the pressure gradients. This reasoning is reinforced by the observation that 
is almost real: the real part describes the "springiness" of the plasma, whereas the small 
imaginary part describes damping. 

We can go further and obtain an expression for the elliptic flow coefficient V2 normalized 
by the eccentricity of the initial state at t = 0. To this end, let us explicitly consider the 
50(3) rotation that rotates the z axis into the y axis. After this rotation, there are two 
more modes involved: 



R-s = {scalar, n = 3, i = 2, m = 2, pn/L = 13, 1.637 - 0.064i, p = 0} 

Ri = {scalar, n = 3, i = 2, m = -2, ph/L = 13, Q ^ 1.637 - 0.064i, p = 0} . 

Consider an expansion 

e = (Too) = eo + ei + ... eo = eS^\t,r) 



;iio) 



;iii) 



where . . . indicates other deformations that we do not include in subsequent calculations. 
The coefficients corresponding to a 90° rotation of (11091) are 

Q3 = Q4 = -\ Iq2 Q2 = -Iq2. (112) 

V o Z 

The eccentricity at time t = is 

^ _ {y"^ - x^) _ (r2sin^6'cos20) ^ /j^3ei''^^sin^6'cos20 . 
~ + (r2sin^6') j-^^ r'^ sir? 9 

The averages are computed with respect to a probability measure proportional to the energy 
density e given in (lllip . evaluated at t = 0. The final expression in (I113P is approximate 
because we used eo in the denominator instead of e = e = eo + ei. 

We will define V2 to be (cos 20) evaluated at = 7r/2 (mid-rapidity), where the average is 
again computed with respect to energy density, but now evaluated in the limit of late times. 
The energy density e is known in closed form, and it turns out that the integrals that define 
the averages can be done. The result is 

— = — Re <^ — sm — ) . (114) 
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Readers wishing to verify (11141) explicitly will be aided by the formula 
e^s = je^ I — -h I sm^ 9 



IT [r4 - 2r2(-l + t2) + (1 + t2)2]4 \^^2 _ _ 



X 



r\-40 - 20int - + 3nH'^) - 2r\A0 + Q4t^ + QH"^ - 8t^ + 3nH^) 

2n2/ An I on„-Oj- oj.2 , oo2j.2^ 



;i + ty{-Ao + 20int - sr + ai^V) 

(115) 

The Q corresponding to the low-lying n = 3 scalar mode is not very sensitive to the choice 
of ph- see table El Plugging in the value for Q quoted in R2, R3, and R4 gives V2/S = 0.37. 
This is at least in the right ballpark for comparison to experiment. For example, in [55], after 
PT^-averaging the elliptic flow coefficient over 0.3 — 2.5 GeV/c and estimating the eccentricity 
using a Glauber model, a value 5/t'2 ~ 3.1 is found. 
Weak points in our analysis include: 

1. The formula (11141) appears highly predictive, but we should remember that there are 
d-wave excitations for higher n as well. Linear perturbation theory by itself cannot say 
whether the minimal n = 3 rf-wave mode is the dominant one. This is one of many 
aspects where a better understanding of collisions of black holes in global AdS^ (for 
example through numerical simulation) would be desirable. 

2. The configuration at t = is nearly spherical for small eccentricity 6, whereas in 
heavy-ion collisions at RHIC, the t = configuration for small eccentricity is Lorentz- 
fiattened into a disk. Again, this is curable in principle if one could numerically solve 
the full Einstein equations for black hole collisions in AdS^. 

3. We have neglected hadronization. This could be fixed simply by appending a Cooper- 
Frye prescription [HH] to the evolution (II 111) of the energy density. 

4.3 Fast modes and the thermalization time 

Aside from the slow modes discussed in sections 14.11 and 14.21 all deformations of the radial 
conformal soliton damp out quite quickly. The time dependence of these modes is dominated 
by the factor e~*^^/^. Referring to (1501) . we see that r ^ 2t near the origin, t = r = 0. So 
for |r| ^ L, the time-dependence of the fast modes is roughly g-sim/L^ ^^^^ conclude 
that the timescale for their decay by a factor of 1/e is 
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To obtain the first estimate in flll6p we used the frequency Qi in the main sequence with 
the least negative imaginary part. To obtain the second approximate equality we used fl57p . 
To obtain the third approximate equality we plugged Tpcak = 300 MeV, which is in the 
upper ranges of temperatures believed to be attained in thermalized states at RHIC. It is 
natural to conjecture an extrapolation of the result flll6p to more general conformations of 
A/" = 4 plasma: l/rg-foid ~ 8.6Tiocai for the time Te-foid that it takes deviations from local 
equilibrium to decay by 1/e. 

The thermalization time Ttherm in RHIC collisions is believed to be in the vicinity of 
0.6 fm/cl^ Larger values tend to cause elliptic flow to be under-predicted by hydrodynami- 
cal simulations. But it is a bit of a puzzle in QCD how such a fast equilibration time results. 
Analysis of perturbative scattering processes seem to indicate rthcrm ~ 2.5fm/c |58l [59] . 
The momentum scales of typical processes decrease quickly during the approach to ther- 
mal equilibrium in RHIC collisions, so the effective 't Hooft coupling quickly increases away 
from the perturbative regime. But it is not clear that the regime of couplings pertinent to 
equilibration is directly comparable to the A — ^ cx) limit of A/" = 4 SYM where supergrav- 
ity gravity calculations become reliable. So we should be more cautious than usual about 
attempting to apply flll6p to RHIC collisions. Nevertheless, it is comforting to see a short 
e-folding time in a strongly coupled plasma. To put it another way: it would have been 
alarming to find Tg-foid > 1 fm/c from strongly coupled A/" = 4 SYM, where if anything one 
should expect faster thermalization than in real-world RHIC collisions. Given the highly 
anisotropic momentum space distribution expected in the early stages of a RHIC collision, 
perhaps it is reasonable to expect that several e-folding times of the relevant thermalization 
/ isotropization processes must elapse before hydrodynamic approximations can be used. If 
so, then perhaps starting from (11160 we should estimate rthcrm ~ 0.3 fm/c for the analo- 
gous processes in A/" = 4 SYM. One should bear in mind that this estimate incorporates 
considerable uncertainties. 

One of the favored thermalization scenarios in RHIC collisions involves plasma instabil- 
ities (see for example the review talk [60] and references therein), which are argued to be 
capable of producing a spatially isotropic form of the local stress tensor in a time comparable 
to 0.6 fm/c [57]. The plasma instabilities owe their existence to an anisotropic distribution 
of momenta and have their wave-vectors primarily along the beamline. Perhaps some of 
the main sequence quasinormal modes describe the late stages of the evolution of analogs of 

^^In [57] the value Ttherm ~ 0.6 fm/c is quoted, and in [5] the range 0.6 fm/c < rtherm < Ifm/c is favored. 
The basis for this range of values appears to be a survey of hydro calculations with various calculational 
methods as well as differing assumptions about the equation of state. 
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these plasma instabilities. 

5 Discussion 

Through the conformal transformation fjlSl) we have mapped a static plasma on S''^ x R to 
a radial flow of thermal matter in Minkowski space which we have termed the "conformal 
soliton." Our subsequent calculations of quasinormal modes in the GAdSBH background 
can be interpreted holographically as slight motions in the static plasma on S^, or as per- 
turbations to the radial flow. A salient feature of these perturbations is that they oscillate 
and damp out exponentially in global time r, which runs over only a flnite range of values 
across Minkowski space. Indeed, the qualitative physics behind our estimate V2/S ~ 0.37 is 
that a d-wave excitation has about enough global time to go through half an oscillation. 

5.1 Linearized hydrodynamics 

Given the separation by a factor of tens to hundreds (depending on the choice of pn / L) be- 
tween the imaginary parts of the low-lying scalar and vector modes and the "main sequence" 
modes, it is natural to regard the former as hydrodynamic and the latter as microscopic. 
With this in mind, we might expect to reproduce the low-lying modes from linearized hy- 
drodynamics. It is easiest to think about linearized hydrodynamics on x R: the general 
equations are 

V'^T.fe = 

fab = (e + P)UaUb + PQab + Tab 

( ~ ~ 2 ~ \ ~ (117) 

Aab = (jab + UaUb , 

and the condition = implies both e = 3p and = 0. As usual, u"' is a timelike 
vector with = —1. The metric cjab is the one specifled in fl32|) . and Va is the associated 
connection. A slight subtlety arises in deflning Tab because of the presence of the Casimir 
energy, which makes a temperature-independent pure trace contribution to the stress tensor. 
We will ignore this contribution and equate Tab with {Tab)suh as appearing in and deflned 
in the discussion preceding (!33l) . so that the unperturbed fluid has stress energy as given in 

dSD. 
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In the linearized approximation, Ua = (—1, Ui) where i runs over the directions and 
the Ui are smalL One also must write a perturbed expression for the pressure, p = Po + Sp, 
where 6p is the same order as the Ui. The resulting equations are 



3dr5p + ApoW'ui = 

2r] 



:il8) 



ApodrUi + di6p - 7] (v^VjUi + V'V%^ + -jdiiV^Uj) = . 

To ease the notational burden, let us set L = 1 for the rest of the discussion. 
An appropriate ansatz to describe the low-lying scalar modes is 

6p = Kie-'^^^ Ui = i^se-^^'S, , (119) 

where S is a scalar harmonic on 5*^ and Sj = — ^(9jS as in (12T!) . Plugging (11191) into (11181) 
leads to 

-SiQKi + ApoksK2 = 
-ksKi + -iipoQ + -^kl -4r]]K2 = 0. 



3 

One may eliminate Ki and K2 to obtain a quadratic equation for Q. For small 77, the roots 
are 

n = ±^^ 'S>i_^±^oif). (121) 

V3 6 Po 

To simplify (11211) . we use po = eo/3 and then combine the relations ([7]) for the total entropy 
S and mass M of the unperturbed GAdSBH solution into the form 

H = ^1. = ^''^ (122) 

Po s M s 1 + pI' ^ ' 

To leading order in large pn, 

n - ±^ - ^(^^ - (123) 

where we have used (11221) and set t]/s = I/Att. Equation (I123p is in good agreement with 
the results of numerical evaluation of low- lying scalar quasinormal frequencies for large pn'- 
see f lTOTD . ffT02D . and figure O 

An analogous treatment can be made for the low-lying vector modes, where the appro- 
priate fluid dynamics ansatz is 

Sp=0 Ui = Kse-'^^Yi . (124) 
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The first equation of (5.2) is satisfied trivially, and the second leads immediately to a linear 
equation for Q 

-4ipofi + r/ (A;^ - 2) = (125) 



whose root is 



n = -i'.^^^^. (126) 

4 Po 



Using (5.6) and rj/s = l/47r, one finds to leading order in large pn that 

n = -^^, (127) 

which is in good agreement with the fit (3.43) to our numerical evaluations of low-lying 
vector quasinormal frequencies. 

A hydrodynamic ansatz cannot be constructed from the spherical harmonics Tij and their 
derivatives because of the conservation and tracelessness properties of Tij. This tallies with 
the absence of low-lying tensor quasinormal frequencies. 



5.2 Boosted conformal solitons 

As we have previously remarked, a natural extension of the linearized quasinormal mode 
calculations presented in this paper would be a finite-element analysis of colliding black 
holes. To gain some intuition for what such an analysis might show, we present the results 
of two small calculations: first, in this section, the holographic "shadow" (T^n) of two black 
holes about to collide; and second, in section 15.31 the geometry in AdS^ of massive test 
particles on near-collision trajectories. In section [5^ we make some heuristic remarks about 
black hole collisions. 

Recall that the conformal soliton describes a spherical shell of conformal matter that 
first implodes and then expands. It is well understood that a Lorentz boost of a spherically 
symmetric shell expanding at the speed of light is still a spherical shell: only the energy 
density is no longer even distributed over the sphere, so that the center of mass moves. 
This does not seem very like a gold nucleus moving at some speed v. But as v 1, most 
of the stress-energy is concentrated near the "beamline," and the boosted radial conformal 
soliton begins to resemble a Lorentz-flattened pancake, with transverse size roughly 2L and 
thickness roughly 2L/7. See figure El 
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A) t = -1.5, V = 0.25 B) t = -1.5, v = 0.75 C) t = -1.5, v = 0.99 



X 



-2-1012 -2-1012 -2-1012 

z z z 

Figure 6: Energy contours of boosted conformal solitons. L has been set to 1. The red 
hne indicates the "beam hne," parallel to the center-of-mass momentum. (A) shows that 
at slower velocities one starts to recover the picture of radial inflow. In (B) and (C) we 
have indicated an additional boosted conformal soliton (unshaded contours) moving in the 
opposite direction, with zero impact parameter. The energy density is azimuthally symmetric 
around the beam line. 

5.3 Test masses in AdS^ 

It is interesting to recall that the conformal soliton depends on the size ph of the black 
hole only through a normalization. So exactly the same contour plots as shown in figure [6] 
would result from considering the holographic shadow of very large or very small black holes 
moving through AdS^, only with a different scale for the energy density. In the small pn 
limit we may replace the black holes with test masses moving on timelike geodesies in AdS^. 
All such geodesies are 5*0(4, 2) transformations of the trivial one, p = 0, representing a test 
mass at rest (and hence an unboosted conformal soliton with very small amplitude). Let 
us trace through the steps to describe the test mass trajectories corresponding to a head-on 
collision with each mass moving at speed v: 

1. In Minkowski space, the trajectory of the center of mass of the right-moving soliton is 
described by the equations 

= vt = = . (128) 

2. Using (1421) and our knowledge of the special case f = 0, we reason that the appropriate 
trajectory of the test mass in AdS^ is described by 

= vX° X^ = X^ = X^ = 0. (129) 
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This trajectory, like any timelike geodesic in AdS^, is the intersection of the hyperboloid 
f l39p with a two-dimensional plane (in this case fll29p ) that passes through the origin. 



3. Using (140|) . one obtains from fll29p the following description of the test mass trajecto- 



ries: 



p/L 



T 

V sin — 



TT 



At r = 0, when the test mass corresponding to the right-moving conformal soliton 
passes through p = 0, its 9 switches from vr to 0. The opposite happens for the left- 
moving test mass. Note that there is only one intersection point of the two geodesies 
in the wedge of AdS^ covered by Poincare coordinates (see figure [TK). 

4. Using (H2l) directly, or (jHj) and f ll30p . one may show that the trajectory of the right- 
moving test mass in the Poincare patch is described by (11280 together with 

y = Vl + (1 - t;2)t2/L2 . (131) 
IJ 

Now consider a non-central collision of boosted conformal solitons. We start by offsetting 
each one in the x direction by ±6/2, where h is the impact parameter. This is trivially 
implemented in the Poincare patch description: the trajectory of the right-moving test mass 
is 

x^ = vt x^ = \ ^2 = ^ = v/l + (l-t'2)tVL2, (132) 

and for the left-moving mass one reverses the signs of v and h. Using fH21) one immediately 
finds the description of the trajectory fll32p in terms of the hyperboloid (l39l) cut by the plane 

X^ = vX'' X' = —{X-' + X') X^ = X^ = -—X\ (133) 

27L AL 

Using (Ho!) one may rewrite fll33p (after some algebra) as 

p/L V sin Y AL ^ Ah cot f 

tanx = -r— ^ tane = -— (134) 



V^l + pV^^ sinxcos6' 6sin6' vLS + b'^/L'^ 

It is straightforward to solve (11340 explicitly for p, x, and 6 in terms of r, but the resulting 
expressions are long and not very enlightening. 
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5.4 Colliding black holes in AdS^, 

Black hole collisions are understood to be highly inelastic, with a large fraction of the total 
energy going into a single final state black hole. The collision of black holes much smaller 
than the AdS^ radius {ph/L <^ 1) could be described approximately as a collision in flat 
space. The four-dimensional case is a classic problem of numerical general relativity, and 
considerable progress has been made on it recently: see for example [nH E2]- To make the 
closest possible analogy with RHIC collisions, one should consider highly relativistic collisions 
of black holes with p// ^ 1 in five dimensions — a significantly harder problem! But unlike in 
the astrophysically relevant case, one does not need to evolve for a long time: global time need 
only run from r ^ —ttL to ttL. The remarkable separation in Minkowski space of conformal 
solitons boosted in opposite directions, illustrated in figure [6], suggests that when the speed 
V is close to 1, even large black holes are in some sense well-separated and non- interacting 
in AdS^ until a global time |r| ^ L, possibly as small as 2L/7. However, it is understood at 
least in the flat space case that a horizon forms around colliding shock waves that encloses a 
large percentage of the energy [63l [Ml ESI [66l [67]. For purposes of describing the subsequent 
"ring-down," the geometry outside the newly formed horizon is more relevant than the 
microscopic interactions inside it because the latter are, by definition, causally inaccessible 
to an outside observer. A thought-provoking idea is that some holographic translation of this 
picture should be relevant to heavy-ion collisions. The possible relation of horizon formation 
from colliding shock waves in AdS^ to RHIC physics has been explored already in [TT f [T8l [T9] . 
Heuristically, horizon formation means that perturbative processes early in the collision are 
washed out, or masked, by subsequent strongly interacting dynamics. For certain high pt 
processes, such as the production of direct photons with px ^ 4GeV/c, this can hardly be 
true. Nevertheless, it seems to us plausible that horizon formation in AdS^ might have to do 
with opacity properties of the QGP. Clearly it would desirable to supplement these highly 
speculative notions with calculations relevant to RHIC observables. 

5.5 Summary 

We have shown how to describe an expanding plasma of finite extent in terms of the holo- 
graphic image of a global AdS^ black hole (GAdSBH). The key step is to use a conformal 
transformation to map a patch of S*^ x R (the boundary of global AdS^) onto R"^'^. Before 
this transformation, the unperturbed black hole is dual to a plasma at rest on S^. After 
the transformation, the unperturbed plasma flow is entirely radial, first contracting and then 
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expanding. In order to describe a plasma whose extent (at the time of maximal compression) 
is much greater than its reciprocal temperature, we must consider black holes whose horizon 
radius pu is much larger than the scale of curvature L of AdS^. 

The bulk of our calculations are, as far as we know, the first full treatment of gravitational 
quasinormal modes of the five-dimensional GAdSBH background. These perturbations are 
dual to slightly anisotropic expansions of the plasma. For large black holes, there is a clean 
separation of modes into hydrodynamic and microscopic modes. The microscopic modes 
decay with an e-folding time that is comparable to the reciprocal of the first Matsubara 
frequency, suggesting a thermalization time Ttherm ~ 0.3fm/c for temperatures and initial 
anisotropics comparable to those at RHIC collisions. Analysis of one of the hydrodynamic 
modes leads to an estimate V2/5 ~ 0.37, where V2 is the elliptic flow coefficient and 5 is the 
eccentricity of the initial state. 

The reader is emphatically warned that applications of our calculations to heavy ion 
collisions rest on assumptions and uncontrolled approximations. Among them: 

1. We have replaced QCD by A/" = 4 super- Yang- Mills theory in the limit of strong 
coupling. 

2. We rely particularly heavily on conformal invariance in passing between S*^ x R and 
Minkowski space. 

3. The unperturbed motion of the plasma is radial rather than longitudinal. 

We have already remarked that point [3] could perhaps be ameliorated by better analyzing 
collisions of black holes or shock waves in AdS^. Our quasinormal mode analysis was entirely 
based on linearizing Einstein's equations around a globally static background, and it does not 
seem likely that QNM's are wholly adequate to describe the evolution of a horizon that forms 
from a high-speed collision. However, the qualitative successes of the simplest case studies 
in section |4] suggest that QNM's will continue to play a role even once more sophisticated 
analyses or full-fledged numerical studies become available. 
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A Spherical harmonics on 



Scalar, vector, and tensor spherical harmonics on are well studied, in part because of 
their utility in describing perturbations of the closed four-dimensional FRW universe |68j . 
Other early treatises include [691 1701 [50] , and a more general account of harmonics on for 
conserved, symmetric, traceless tensors of rank r was given in [7T]. Mostly we follow |50j . 
We follow standard notation for scalar spherical harmonics on S"^: 



YeUOA) = W^|:p^^^(cose)e^-^ (135) 



where the associated Legendre polynomials are 



pn^) = - -'r"i^(^' ~ 1)', (136) 



so that, for example, Yoo = l/v4vr and 

Fn = -J— sin^e*^ Yio = \-cose ^ _i = J— sin ^ e"*^ . (137) 
V Svr V TT ' V Svr 

Also, it is convenient to introduce functions 

*,..W = ^^^, (138) 



d{cosx) 



which are non-zero for i < n. 
Spherical harmonics 



U(nim) = <l>ni{x)YUO, 4>) a^^im = ^T"^^^^ (139) 

^JCinim 2 [n-i-l)\ 

obey (IE]) with kg = — 1, where |m| < i < n. The E>{nim) for fixed n form n^-dimensional 
representations of S'0(4) = SU{2)l x SU{2)ji, namely spin-(ra — l)/2 for both SU{2) factors. 
The Snim are orthonormal: 

/ sin^ X sin dx dO dcf) 'S{n£m) S* {n'i'm') = 5nn'^u'^mm' ■ (140) 
J53 
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The first few E>{nim) are 



\/2 

§(100) = 1/7tV2 §(200) = — cos X 



1 /2 1 

§(211) = sin X sin ^e''^ §(210) = — sin^cos^ §(2, 1, -1) = - sin x sin ^ e"*'^ . 

71 71 71 

(141) 

Tfie group tlieory for harmonics on of conserved, symmetric, traceless tensors of rank 
r is simple to state. (Vector spherical harmonics are the special case r = 1, and tensor 
spherical harmonics are the case r = 2.) Let ul denote the n-dimensional representation of 
SU{2)l (i.e. spin {n — l)/2), and similarly for SU{2)r, so that the E>{nim) for fixed n fill out 
the representation ni ® tir. Then, for fixed n, defined so that the eigenvalue of the laplacian 
is —k^ = —n^ + r + 1, with integer n > r, the representation content of the harmonics with 
that eigenvalue is 

{n + r)L ^ {n - r)R © (n - r)/, (g) (n + r)^ , (142) 

for a total dimension of 2(?7,^— r^). But there is a Z2 parity symmetry under which SU (2) l and 
SU{2)r are interchanged, and the representation fll42|] may be decomposed into symmetric 
and anti-symmetric product pieces, each with dimension ra^ — r^. 

Explicit expressions for vector spherical harmonics were given in [50] and are reproduced 
here. For the even modes (corresponding to the symmetric part of the product fll42p ). 



Vf = ^^,Vr Vr = \l ^-^^ sin-^ X $n.(x) 
t{i + 1) y n^anim 



( vf \ 

Vfde 



(143) 




where anim is defined as in fll39p and Q < i < n. The eigenvalue for these modes is 
ky = 'n? — 2. The components of the vector are in the Xi <^ directions, and the index 
is down, meaning that Vidy^ is a one-form. Thus, for example, = means that 
Yidy'' = dx- The V^™'^(n£m) are orthonormal: 

sin^ X sin 6 dx dO d(j) V"™"^'^' [nim) V"™'''* {n'i'm') = <5„„'<5«'<5w , (144) 

where the index j on the first factor is raised using the standard metric on the unit 5''^. The 
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first few V°™'^(n£m) are 



yeven(^211) 



IT 



V2 



sin 6 



\ 



sin X cos X cos 9 
\i sin X cos x sin 6*^ 



Vf"°(210) = - 



IT 



( cose \ 

— sin X cos X sin 6* 




ycvcn^2, 1,-1) 



-i(f> 



TX 



V2 



sin0 



(145) 



sin X cos X cos 6 
\—i sin X cos x sin 



Tlie odd vector modes are 



Vf'^(n£m) 



sin X ^ / s 1 



l)anem 



sin 9 



( \ 



(146) 



where we have altered the overall normalization from the one stated in [30]; which appears 
to be slightly in error. Again, < £ < ra, and the eigenvalue is ky = — 2. The Yj'^'^{n£m) 
are orthonormal: a formula precisely analogous to (I144p applies. And they are orthogonal 
to the Vf™(n£m). The first few Yf'^inim) are 



/ 



j0 

Vf ^(211) = — ^ 



I sm X 



( 



¥^^(210) 



^— sin^ X sin ^ cos d ^ 



1 

TT 








\ 



ysin^ X s\\? b j 



( 



(147) 



¥^^^(2,1,-1) 



TT 
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% sm X 



ysin^ X sin Q cos 6* y 



It is easy to verify that there are r? — \ even and r? — \ odd vector harmonics at level n > 2, 
as the representation theory (11421) demands. 

Even tensor harmonics, again following [5U], take the form 



T^™°(n£m) 



Tf{de - cot 



\T^% Tf{de - cot e)d^ Tf{dl + sin 6 cos Ode) + Tf sin^ 6 j 



(148) 
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(149) 



where 

Tr=^|^(5x + 3cotx)Tr 
^3"- + 2 cot x)T.-T-] 

= ^ + 1)T3"^ - sin' xT^rT 

and n > I > 1. The eigenvalue for these modes is k"^ = — 3. The T^™'^(n£m) are 
orthonormah 

/ sin' X sin ^ rfx dO d(j) r™'^'^^' (n£m) T^™'''* {n'i'm') = 5„n'5u'5mm' , (150) 

The = 3 even tensor spherical harmonics, ¥^^^^(3, 2, m), fill out a 5-dimensional represen- 
tation of 5*0(4). The explicit form of these tensors is too complicated to be worth displaying 
explicitly here. 

Odd tensor spherical harmonics take the form 

T°/<^(n£m) = 

f -TfcscOd^ Tfsinede ^ 

-Tf CSC e -2Tf CSC e{de - cot e)d^ Tf{^\x^ ^ 9| - cos Qde- esc Q 3%) 

TfMde Tf{?>inedl-co^ede-csc6dl) 2T^^ sin 6 {de - cot 6)8^ j 



where 



(151) 



^."'= (,-i)Wi) ^''-'"^^^"^- 

Again n > £ > 1, and = n"^ — 3. We have altered the overall normalization on T°!^'^{nim) 
from the expressions given in [50], which appear to be slightly in error. The T°^^(n£m) 
are orthonormal, and they are orthogonal to the T°™'^(n£m). The expressions for even the 
low-lying T°^'^{ni'm) are too complicated to reproduce here. As with the vector modes, it is 
easy to verify that there are n' — 4 even and n' — 4 odd tensor harmonics at level n > 3, as 
the representation theory (11421) demands. 
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